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The fabrication of more and more miniaturized electronic and photonic devices rehes on new, 
ingenious methods for the fabrication of spatiahy controhed nanostructures. Examples are elec- 
tronic devices based on semiconducting nanowires and photonic devices based on chains of metallic 
nanoclusters that guide light by coupled surface plasmons. In this contribution a template-directed 
ion beam synthesis of nanowires and regular nanocluster chains will be presented. As templates, 
V-grooves etched in (OOl)Si and subsequently oxidized are used. High fluence Ge^ implantation is 
carried out into the Si02 layer at 70 keV. Thereby, the implanted Ge enriches themselves in the V- 
groove bottom to a critical amount, which may result in nanowire formation by nucleation, growth 
and coalescence during subsequent thermal treatment. TEM investigations indicate the formation 
of a nanowire buried in the Si02 at the V-groove bottom. 

Kinetic lattice Monte Carlo simulations of the nanowire formation process were performed in 
order to understand the phase separation mechanism and results are compared to TEM images. 
Furthermore, it is shown that even ideal nanowires show a instability and form during long-lasting 
anneaUng equal-spaced and equal-sized nanoclusters ("nanocluster chains") by self-organization. 



I. INTRODUCTION 

A wide field for novel applications opens at the bot- 
tom of the length scale accessible by modern semicon- 
ductor preparation techniques. Synthesized nanos- 
tructures show interesting physical properties being 
completely different from the bulk material. Be- 
sides other applications, this could help to miniatur- 
ize modern optical devices. Thus, having metallic 
nanoclusters arranged in a regular chain, light can 
be transported by the coherent coupling of surface- 
plasmon-polaritons Along such a chain, light can 
be guided around corners well below the classical 
diffraction limit, allowing the integration of optical 
devices on the chip-level (Plasmonics) [|| . 

However, the precise synthesis of the nanostruc- 
tures needed for such applications remains a great 
challenge to current research. This becomes even 
more important since classical structuring fails at this 
level. Here, self-organization processes might be able 
to assemble structures with nanometer dimensions. 
In particular, the ion beam synthesis (IBS) became 
successful for fabrication of buried nanostructures. 
Quantum dots were formed under various conditions 
H, Using very high ion fluences, buried layers were 
synthesized [^. 

In order to obtain spatially well-defined nanostruc- 
tures, taming of self-organization of pattern forma- 
tion is an attractive approach. In the present contri- 
bution, this is achieved by oxide covered Si V-grooves 
as templates, which guide the self-assembly of buried 
nanowires from implanted Ge atoms during thermal 
treatment. This scheme is depicted in Fig. |l|. By 
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FIG. 1: Scheme of the ion beam synthesis of a buried Ge 
nanowire at the bottom of an oxidized silicon V-groove. 
(a) Ge^ ions are implanted at high fluence into the V- 
groove. Atoms are sputtered from and are redeposited 
onto the sidewalls leading to a mass transport towards 
the V-groove bottom. During annealing (b) in an inert 
ambient a nanowire may form at the V-groove bottom. 



a similar approach Si nanowires were recently fabri- 
cated by 0+ implantation into Si V-grooves |^]. How- 
ever, the method of ref. is restricted to Si wires 
and is hard to control. Here, Ge was implanted as 
it is a well-known material for nanocluster synthe- 
sis. However, also other ions can be easily applied, 
e.g. Au for plasmonic devices. Different stages of the 
fabrication process were monitored by cross sectional 
Transmission Electron Microscopy (X-TEM) investi- 
gations. 

Theoretical studies are presented to understand de- 
tails of the phase separation mechanism. On the 
atomic scale, within kinetic 3D lattice Monte Carlo 
simulations (KLMC), the diffusion of Ge atoms, their 
nucleation to clusters, cluster ripening and coales- 
cence are tracked. For long lasting thermal treatment 
such simulations also describe the pearling instabil- 
ity of wires. A nanowire tends to decay into droplets 
driven by the minimization of interfacial energy and 
highly regular chains of equal-sized and equal-spaced 
nanoclusters form. 
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FIG. 2: XTEM images of (a) an as-oxidized and (b) an 
as-implanted sample (70 keV, 1 x 10^'^ Ge+cm"^). The 
Ge content in the Si02 appears as dark shadow in the 
TEM image (b). 



II. EXPERIMENTS 

V-grooves aligned along the [110] direction on 
(001) Si wafers were prepared by photolithographic 
masking and anisotropic etching with 30% KOH at 
80 °C. This self-adjusting process results in atom- 
ically smooth (111) crystal facets forming the side- 
walls of the V-grooves with an angle of 54.7 ° to the 
plane wafer surface. V-grooves of 4 \xm width and 
5 mm length were studied. They were oxidized in 
dry O2 at 1000 °C, leading to 200 nm thick Si02 on 
the (OOl)Si surface and 220 nm thick Si02 on the 
(111) sidewalls of the V-groove. A sharp surface de- 
pression forms during oxidation at the V-groove bot- 
tom, which can be seen in X-TEM images (Fig. ^ a). 
The Si02 covered V-grooves were implanted at room 
temperature with 1 x 10^^ Ge+cm~^ at 70 keV (flux 
density 3 x 10^^ cm~^s~^). The content of Ge in the 
Si02 appears thereby as dark shadow in the TEM im- 
age (Fig. Ijb). During implantation, Ge is enriched in 
the V-groove bottom compared to the sidewalls due 
to (i) the reduced ion density per surface area at the 
sidewalls caused by the inclination and (ii) due to ion 
erosion of the sidewalls. Sputtered atoms including 
Ge are redeposited at the opposite sidewall and refill 
the sharp kink at the as-oxidized V-groove bottom. 
The V-groove surface evolution under ion bombard- 
ment was recently modeled in a continuum framework 
in nice agreement with V-grooves shapes observed 
in TEM images also proving a Ge accumulation of 
30% in the V-groove bottom. 

Subsequent furnace annealing in N2 ambient was 
performed at 950 °C for 10 minutes. Thereby, the 
small amount of implanted Ge is extremely sensitive 
to moisture, which may result in the oxidation of 
Ge. Therefore, (i) the standard wet cleaning step 
in H2O2/H2SO4 prior to any thermal processing was 
skipped and (ii) the samples were transferred into the 
cold furnace and then purged in cold N2 in order to 
reduce the moisture intake. Fig. ^ shows an XTEM 
image of an annealed sample. Buried Ge clusters are 
present in the sidewalls and the bottom. The dark 
dot, visible in the V-groove bottom, is a cross-section 




FIG. 3: XTEM image of an annealed sample. The large 
dark spot in the V-groove bottom appears to be a cross 
section of a Ge nanowire with 35 nm diameter. 




FIG. 4: High resolution TEM image of a Ge cluster in an 
annealed sample showing lattice fringes, which proves the 
crystallinity of the precipitates. 



of the expected nanowire with a diameter of 35 nm. 
All Ge cluster are crystalline, which is confirmed by 
a high resolution TEM image. Fig. ||, of a Ge cluster 
showing lattice fringes. 

The continuity of the wire over long distances re- 
mains to be proven. Several XTEM images taken 
from different cross-sections show always the large 
black dot with identical diameter and position in the 
V-groove bottom, which is an indication for a contin- 
uous nanowire. 



III. KINETIC MONTE-CARLO 
SIMULATIONS 

KLMC simulations were performed in order to 
understand the wire formation from dispensed Ge 
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atoms. The interactions between impurity atoms 
(Ge) were simulated by a Cellular Automaton ap- 
proach. Gauging the bond strength Eb of Ge atoms, 
their effective nearest-neighbor interactions within 
the host matrix, the Si02, must be considered. Here, 
the most simple implementation is used, i.e. the in- 
teraction is assumed to be configuration independent. 
The binding energy uEb of a Ge atoms is just pro- 
portional to the number n of nearest-neighbor sites 
occupied by Ge and can be described by the Ising 
model. Additionally, the diffusion of Ge monomers 
in the host matrix is thermally activated with an ac- 
tivation energy Ea- 

Then, the transition from the initially out-of- 
equilibrium state towards thermodynamical equilib- 
rium is described by the importance sampling of con- 
figurations according the Metropolis algorithm ||] . A 
sequence of states is produced representing the time 
evolution, where each state is generated from its pre- 
decessor by a jump of a randomly chosen Ge atom 
to an arbitrary empty nearest-neighbor site with the 
transition probability 



exp |— (n- — ^ 



Ea 



I , Uf < Hi. 



During one Monte Carlo step (MCS) each atom per- 
forms statistically one attempt to jump from its ini- 
tial site i to a final site /. Thereby, the activation 
energy for diffusion Ea is assumed to be the same for 
each possible transition, i.e. for interface and bulk 
diffusion. Thus, it is possible to renormalize the tran- 
sition probability Pif such that each diffusional jump 
is successful {Ea 0), which only scales the time of 
a MCS thereby improving the simulation speed 
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exp{-(ni-n/)0T}, 



Uf < Hi. 



This approach has been proven useful for the descrip- 
tion of precipitation after ion implantation and 



the evolution of fee surfaces under ion erosion |10 . 

For the KLMC simulations of the nanowire forma- 
tion a simplified initial distribution of Ge atoms was 
assumed. A cylindrical profile with a Gaussian cross 
section aligned to the [100] direction of the simulation 
lattice was used (Fig. || a) with a peak concentration 
of 31% and4 nm standard deviation. 95% of the im- 
planted atoms lay in a cylinder of 16{nm diameter, 
which resembles the Ge-rich region at the V-groove 
bottom, see Fig. |b). 

During annealing, Ge nanocluster nucleate from 
the supersaturation and grow in time. They mostly 
have nucleated at the position of highest supersat- 
uration, which is the profile center. As they grow, 
they start to touch each other and coalesce to ex- 
tended structures. A continuous wire forms if this 
coalescence proceeds further, while the surrounding 
nanoclusters evaporate and Ge condenses onto the 
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FIG. 5: Kinetic lattice Monte Carlo simulation of the 
nanowire formation from an initially cylindrical Gaussian 
distribution of impurity atoms with a peak concentration 
of 31% and 4 nm standard deviation. Shown are (a) the 
initial, (b) an intermediate (102 kMCS), and (c) the final 
(490 kMCS) state. The final wire diameter is approxi- 
mately 5 nm. 



wire-like structure, see Fig. ||b,c). The complete co- 
alescence occurs only, if the initial profile is not too 
broad and the Ge concentration not too low, i.e. if 
peak concentration of the initial Ge profile is larger 
than 25% (at the given profile width). Otherwise, the 
clusters remain isolated and randomly distributed in 
size and space. 

Comparing experiment and KLMC simulations, it 
is obvious that a buried nanowire could form at the 
V-groove bottom at the chosen Ge+ ion fiuence. As 
seen from Fig. ^, the annealing has been interrupted 
before all surrounding nanoclusters have disappeared, 
which compares to Fig. |5|b). 

It should be noted that the KLMC simulation uses 
a lattice. Thus, a wire formed by coalescence will 
be single crystalline, whereas an experimentally fab- 
ricated wire is expected to be polycrystalline. There- 
fore, the infiuence of grain boundaries is not taken 
into account in the present study. 



IV. PEARLING INSTABILITY 

For long-lasting annealing, theory predicts that a 
continuous wire decays. Thus, a wire is an unstable 
configuration. Its decay is driven by the minimization 
of the interfacial energy. Harmonic perturbations of 
the wire radius lower its interfacial energy if their 
wavelength is larger than the circumference of the 
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FIG. 6: Kinetic lattice Monte Carlo simulation of the pearling instability of a nanowire with 8.6 nm diameter. The 
wire decay proceeds via interface diffusion and a chain of nanoclusters with an average diameter of 16.4 nm is formed. 



wire, A > 27TR n. T his result has been first found by 
Lord Rayleigh [|ll|] in 1879 for inviscid liquids. Linear 
stability analysis also revealed a fastest growing mode 
of radius fluctuations during the initial stage of 
the fragmentation. Thus, perturbations with wave- 
length Am grow at the fastest rate in time and are ex- 
pected to dominate fragmentation. The value of A^ 
depends only on the specific mass-transport mech- 
anism and not on material parameters or tempera- 
ture. For liquid "wires", which changes their shape 
by hydro dynamic flow, a value of A^ = 9.02 • i?o was 
predicted [|ll[ , which is also valid for internal volume 
diffusion as shown by Nichols and Mullins For 
surface diffusion, it holds 

Xm = 7227ri?o = 8.89 • Rq, 

while for external volume diffusion the value a 12.96 • 
Rq was obtained for A^ [0. With the result of the 
stability analysis an average droplet diameter dduster 
can be given, if one assumes that the volume between 
two consecutive radii minima forms a droplet. For 
surface diffusion this gives 

dduster = \/ 6XmRo = 3.78 - Rq- 

In order to check the stability analysis for large ampli- 
tudes of diameter fluctuations, Monte Carlo simula- 
tions of a wire decay driven by surface diffusion were 
performed. Volume diffusion is explicitly prohibited 
within the simulation, i.e. atoms are not allowed to 
detach from the wire surface and vacancies can not 
be generated within the wire. Snap shots of the wire 
evolution can be seen in Fig. ^. The initial wire, 
Rq = 4.3 nm, shows after a long incubation time sta- 
tistical fluctuations, which grow faster and faster with 
increasing amplitude. They finally lead to the disin- 
tegration of the wire after 5600 kMCS into droplets. 
The average diameter of these droplets is 16.4 nm. 



which is in good agreement with the prediction of the 
linear stability analysis, dduster = 3.78 ■ Rq = 16 nm. 
Note that the simulation was carried out at high 
temperature. The ratio between the effective bond 
strength E^ond and the thermal energy, ksT^ was set 
to 0.7. An absolute value of the temperature T can 
not be given since the value of effective Ge-Ge bond 
strength in Si02, Ei^ondi is unknown. It may be es- 
timated, however, within a good approximation from 
the solubility of Ge in Si02. 

The stability analysis as well as the Monte Carlo 
simulations neglect the influence of grain boundaries 
of a polycrystalline wire. Grain boundary diffusion 
may provide another source of instability, if annealing 
is carried out below the melting point. Then, the size 
of the droplets may differ from the result given above. 

V. SUMMARY 

A novel fabrication method of buried nanowires 
by ion implantation has been demonstrated. Ge+ 
were implanted into oxide covered Si V-grooves lead- 
ing to Ge enrichment in the oxide at the V-groove 
bottom. There, during annealing, the Ge forms a 
nanowire extending along the V-groove. Other than 
a similar method published recently j^], this method 
should work for any implanted ion species immisci- 
ble in Si02. In particular, ion beam synthesized Au 
wires, which are subsequently disintegrated in regular 
nanocluster chains, are of interest for the fabrication 
of plasmonic structures [|[. 

The phase separation mechanisms were studied by 
means of kinetic 3D lattice Monte Carlo simulations. 
As shown, a cylindrical profile with Gaussian cross 
section may coalesce to a continuous nanowire of only 
a few nanometers diameter. For long-lasting thermal 
treatment atomistic simulations revealed the pearling 
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instability. A wire decays into equal-sized and equal- 
spaced droplets. An average droplet diameter of ap- 
proximately twice the initial diameter of the wire, as 
it is predicted by linear stability analysis for a wire 
decay by interface diffusion, could be confirmed by 
Monte Carlo simulations. 

Thus, nanocluster chains having a good periodic- 
ity over long distances can be fabricated by the frag- 
mentation of nanowires using the pearling instability. 
Thereby, the droplet diameter is independent from 
material parameters and is controlled only by the di- 



ameter of the initial wire. 
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